Optimal photon-pair single mode coupling in narrow-band spontaneous parametric 

down-conversion with arbitrary pump profile 
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A theoretical study of the performance of single-mode coupled spontaneous parametric down- 
conversion sources is proposed, which only requires very few assumptions of practical interest : 
narrow-bandwidth and quasi- degenerate collinear generation. Other assumptions like pump beam 
spatial and temporal envelopes, target single-mode profile and size, and non-linear susceptibility 
distribution, are only taken into account in the final step of the computation, thus making the theory 
general and flexible. Figures of merit for performance include absolute collected brightness, pair 
collection efficiency and heralding ratio. The optimization of these values is investigated through 
functions that only depend on dimensionless parameters, allowing for deducing from the results 
the best experimental configuration for a whole range of design choices (e.g. crystal length, pump 
power). A particular application of the theory is validated by an experimental optimization obtained 
under compatible assumptions. A comparison with other works and proposals for numerically 
implementing the theory in its most generality are also provided. 



I. INTRODUCTION 

Sources of entangled photons have applications ranging 
from fundamental tests of quantum mechanics [1, 2] to 
quantum information and communications [3, 4]. Entan- 
glement based quantum key distribution protocols have 
been reported [5] as an alternative to those based on sin- 
gle photon sources and quantum repeaters have also been 
shown to require entanglement as a primary resource 
[6, 7]. 

Spontaneous parametric down- conversion (SPDC) re- 
mains the simplest way to generate entangled photons, 
enabling telecom wavelength generation which cannot be 
easily achieved using atomic cascades. A large choice of 
configurations is available according to the emission ge- 
ometry (colinear [8] or not [9, 10]), the phase matching 
type (I [11] or II [9]) and the signal and idler frequencies 
(degenerate [12] or not [13]). The non-colinear emission 
geometry has some advantages, providing direct separa- 
tion of signal and idler photons but colinear emission in 
periodically poled crystals recently became increasingly 
popular [14-16], because of the high brightness obtained 
with longer crystals and higher nonlinear susceptibilities. 

Filtering is one of the most important aspects of pho- 
ton pair engineering. Spatial filtering and more specifi- 
cally coupling the photon pairs to a single spatial mode, 
like an optical fiber, is the first requirement for long- 
distance communications. The coupling efficiency is 
then a key parameter, knowing that losses directly limit 
the maximum possible visibility of the source [17, 18]. 
Narrow-band spectral filtering can be necessary to en- 
sure spectral indistinguishibility but also to obtain effi- 
cient coupling to a quantum memory [19-22]. It is also 
beneficial for long-distance communications of polariza- 
tion states as it reduces polarization mode dispersion. 
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Since the pioneering works of SPDC concerning the 
quantum fluctuation and noise in parametric processes 
[23] and those dealing with probability of coincidences 
in the emission of signal and idler photons [24-26] an 
important theoretical effort has been developed in order 
to optimize various kind of sources [27, 28], some of them 
dedicated to strongly focused beams [29]. 

After studies of the collection of SPDC through aper- 
tures [30, 31], numerous authors have reported studies of 
the coupling of SPDC into a single spatial mode like that 
of a fiber [32-38], with various assumptions and experi- 
mental methods of validation. 

The optimal focusing of a continuous wave monochro- 
matic pump has been investigated [35] in order to calcu- 
late the maximum coupling efficiency of photon pairs into 
single- mode fibers but the pump diffraction was neglected 
and the phase mismatch was not taken into account in 
the optimization process. The same approximations have 
been used in [36] to calculate the absolute emission rates 
of SPDC into single modes, where the crystal is more- 
over assumed to be thin. The absolute brightness is also 
calculated in [37] but with a continuous wave pump, em- 
phasizing the parallel between SPDC and the classical 
treatment used for second harmonic generation. Very re- 
cently, R. S. Bennink studied the case of Gaussian beam 
profiles, directly writing the nonlinear interaction hamil- 
tonien with a gaussian beam profile for both the pump 
and the SPDC generated photon pairs [38]. In this case, 
neither the absolute brightness of the source nor the ab- 
solute collection efficiency of an optical fiber can be eval- 
uated. Moreover, no spectral filtering is taken into ac- 
count. 

To our knowledge, a comprehensive study of narrow- 
band SPDC pumped by a diffracting beam of arbi- 
trary pulse envelope and duration (from short pulses to 
continuous- wave lasers) and arbitrary transverse profile , 
covering a full range of focusing (from tight focusing to 
parallel beams) and arbitrary filter shapes, has not been 
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FIG. 1. General configuration considered in the theoretical 
framework: a nonlinear crystal is pumped by a beam of ar- 
bitrary temporal and spatial profiles. The photon pairs pro- 
duced through SPDC are spectrally filtered and coupled to 
a single spatial mode, before being split towards channels A 
and B. 



addressed. Such a generalization can, however, be useful 
in view of optimizing future SPDC sources that might 
require specific spectral and spatial characteritics to in- 
terface with particular quantum information systems. 

In this paper, we propose a theoretical framework that 
is general enough to allow the investigation of these char- 
acteristics and their effects on the source performance 
(e.g. brightness, coupling efficiency) under assumptions 
suited to the desirable characteritics mentioned above 
: narrow-bandwidth for quantum information applica- 
tions and colinear emission in long crystals with a given 
non-linear susceptibility distribution. Using a dedicated 
source, we also experimentally validate the theoretical 
predictions for particular assumptions. 

Section II describes the general formalism of the ad- 
dressed problem and derives the wave- function of the cre- 
ated photon pairs, taking into account spatial and spec- 
tral filtering. We also define various figures of merit that 
are useful to characterize a source performance. The nu- 
merical application of this general theoretical approach 
to the case of a narrow-band fibered source pumped by 
a Gaussian beam is developed in section III. The experi- 
mental setup is described in section IV and the measured 
performances are compared with the theoretical previ- 
sions. A comparison of our results with other works, as 
well as guidelines for generalizing the numerical calcula- 
tions are proposed in section V. 



II. THEORETICAL FRAMEWORK 



A. The down-converted two-photon state 

The general configuration considered hereafter is de- 
picted in Figure 1. The non-linear interaction takes place 

_(2) 

in a crystal with a second-order susceptibility % (r), 
pumped by a classical field of positive-frequency com- 
plex amplitude . The spontaneous parametric down- 
conversion process is described by a time-dependent 
Hamiltonian ii{t). The signal and idler fields, and 
E^ , have the following plane- wave decomposition in the 



quantization volume V : 

Eg(r,t) = ^^„a, s ^(4, l ) e i(kf -- r - Wc -* ) (1) 

^ s , i 

where the field amplitude A(£) for the mode i associated 
with annihilation operator hp is 



A(£) 



2e Vn 2 (k e )' 



(2) 



Frequency and polarization are respectively denoted by 
uj£ and §£, wavevectors in the crystal by k^. A notation 
will be used when the wavevector is evaluated in another 
medium. 

In the interaction picture, the Hamiltonian is given by 
[39, 40] 

U(t) = 6 |d 3 rE«(r,t) ^ (2) (r)Ei-)(r,t)E^(r,t) 

+ H.c. (3) 

where is the Hermitian conjugate (H.c.) of E^. 

The pump field is assumed to be a linearly s^-polarized 
classical paraxial beam propagating along the z axis with 
L 2 -normalized temporal and spatial envelopes T p and S p . 

Using transverse coordinates p = ( y ) so that r = ( p z ) 



E< + >(p,z,t) = e x C p T p (t- -)s p {p,z)e l 



■ (4) 



Here k Po is the pump wavevector modulus at its central 
frequency uo Po and v p is the group velocity in the crystal. 
We neglect its spectral dispersion, therefore v p = c/n p 
where n p is the refractive index in the crystal at cu Po . 



The constant C p = y (o p /(2eon p v p ) depends on the pump 

pulse energy § v . 

When the down-conversion efficiency is small, the gen- 
erated two-photon state is the first-order approximate 
solution to the Schrodinger equation. Discarding the 
zeroth-order term 



\m) = ^ / dt"H(t')\0). 



(5) 



Let us separate the spatial dependence R(r) of the non- 
linear susceptibility 



£xX = x(£s,£i) ■ R(r). 



(6) 



After the pump has propagated through the whole crystal 
length, (t — >> oo), the state is time-independent in the 
interaction picture. Identifying the integral in Eq. (5) as 
a Fourier transform of = E^ • e x with respect to its 
temporal variable and according to the Fourier transform 
conventions of Appendix A, the state is given by 



: ii(r)£?«(r,^ a +^JaJ a aJjO). (7) 
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This equation can be understood as the sum of the con- 
tributions of local interactions to the delocalized down- 
converted two-photon state The integral has the 
form of a spatial Fourier transform, which leads to 

IV>> = X>&' £< ) a l a /J°>> (8) 

£ s 

where 

x R*EM(k ia +k ii ,u)i a +ui i ), (9) 

R being the Fourier transform of R. Using the fact that, 
for a diffracting beam in the paraxial approximation, the 
Fourier transform of the transverse spatial envelope has 
the following z dependence 

S p (K,z) = S p (K,0)e~ i ^ z , (10) 

where k is the transverse wavevector, we have, according 
to the Fourier transform conventions of Appendix A : 

££+>(«, k x ,u>) = C p f p (to - to P0 )S p (n, 0) 

x2n5(k z -k P0 -^^ + ^). (11) 

Up ^^po 

We consider a crystal of length L centered on z = 
whose transverse dimensions are large compared to the 
pump beam profile. If its nonlinear susceptiblity distri- 
bution is transversally invariant, it can be expressed as a 
one-dimensional Fourier series 

R(z) = rect(>/L) x ^ i? m e~ i27rmf , (12) 

m 

where rect(z) is the rectangular function, with rect(0) = 
1. In the reciprocal space, 

&( k z) = ^Rm sine (k z + 27rm/L) . (13) 

m 

Then 7 o(4,^) = E m 7o m (4,A) with 

J 0m (£sA) = Tf.RrnX(ts^i)A(£ s )A(£i) 

x CpT p (uj£ s ^uj£ i - Up )Sp(n£ s +k^ 4 ,0) 

x Lsinc ^Aif m (k^ s +ki.,cj £a , (14) 

where 

AK m (k,to) = k z - k P0 - + ^t +m ^. (15) 

Vp ^kpo Li 

So far, we have explicited in Eq.(8) and (9) the un- 
normalized state |^), of a photon-pair down- converted 
during a pump pulse. Its squared modulus (V^) is the 
probability for such a down- conversion to effectively oc- 
cur. The effect of spectral and spatial filtering on the 
photon-pair state is developped in the following section. 



B. Spectral and spatial filtering 

The calculation of the effect of filtering on one-photon 
states is detailed in Appendix B. Extrapolation to two- 
photon states is straightforward and only results are 
given here. Filters give rise to losses and two photon 
states are consequently transformed into mixed states, 
i.e. non coherent superpositions of two-photon, one- 
photon and zero-photon terms. In the case of coinci- 
dence counting, the two-photon state component is post- 
selected and we will discard zero- and one-photon terms. 

Purely spectral filters are described by a function 
^(ijo — uoy) defining their amplitude transmission with 
a maximum normalized at 1 at their central frequency 
cjp- More generally, a spectral filter can be sensitive to 
the direction of the wavevector (like prisms, gratings or 
Fabry- Perot etalons), and it may be necessary to take a 
K-dependence into account : F(uo — ljf,k). 

Spatial filtering is modelled as the coupling of the 
down-converted field to a single spatial mode defined, in a 
medium of refractive index n', by a frequency dependent 
function Wj o(r). The filter location z = zq defines the 
transverse plane where the coupling is considered to take 
place. This coupling induces a state projection on the 
target mode, and the state of the transmitted photons is 
associated with a creation operator 6^ . 

Following Eq. (8), the spectrally-filtered, free-space 
two-photon state is given by 

where 

(2) 

T T ) (£ s ,£i) =J r (u£ 3 -UaoiKle^iuii -Ui ,Kli), (17) 

and uj Sq and are the signal and idler frequencies se- 
lected by the filters, which can be equal for degenerate 
down-conversion using a unique filter. The squared mod- 
ulus of this state, (^ol^o), is the probability Po of gen- 
erating a photon pair during a pump pulse in the filter 
bandwidth. 

The spectrally and spatially filtered state is given by 

£ s ,£i 

(18) 

where 

x e i(fc «.'- + *W*°, (19) 

(9* is the complex conjugate of the transverse Fourier 
transform of 0^,0 > K £ = k Zi in'(oj£)/n(oJi) is the longitu- 
dinal wavevector evaluated in medium of refractive index 
n' , and S is the transverse section of the quantization 
volume. 
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Due to spectral and spatial filtering, has become 
|^2). The squared modulus of this state, (V^IV^), is the 
probability P2 that one pulse has generated a photon pair 
and that this pair has been transmitted by the spectral 
and spatial filters. Its numerical calculation is performed 
in subsection II D, using assumptions that will allow us to 
separate the frequency dependence from the wavevector 
dependence in Eq. (15) as shown in the next subsection. 



C. Assumptions 

A few assumptions are quite natural when aiming at 
applications in quantum information. Moreover, they 
enable further analytical development as well as faster 
numerical calculation. The consequent approximations, 
mostly used to simplify the expressions of AK in (19), 
aim at eventually decoupling the frequency dependence 
(cj) from the angular dependence (n) in the wavevectors 
k. 



1. Collinear collection 

When using long crystals, as in most recent and ef- 
ficient devices [41, 42], a collinear configuration is re- 
quired. In such a case, the paraxial approximation, as 
already used for the pump beam, can be applied to the 
down-converted photons on the z-axis, giving the follow- 
ing expression for the longitudinal component of the sig- 
nal wavevector 



&Zs — ^Sq H~ 



1 \*s\ 2 

2 IkJ ' 



(20) 



where k So is the wavenumber of collinearly emitted signal 
photons at the central frequency of the filter u; So , and 
v s = c/n s with n s the refractive index in the crystal at 
the signal frequency. A similar expression applies for k zi . 



2. Narrow bandwidth 



When filtering limits the source bandwidth to less than 
a few nanometers (Acjp <C ^ Po ), we have uo s w uo SQ and 
&i ~ 0Ji in which case, according to Eq. (20), the phase 
mismatch AK depends only on K s1 Ki. The fields ampli- 
tudes and nonlinear susceptibility, respectively defined in 
Eqs. (2) and (6), can then be considered constant. 

A narrow bandwidth is a desirable characteristic for 
long-distance communications, as it reduces dispersion in 
optical fibers and is also required for projects of repeaters 
based on quantum memories [3, 20]. In such memories 
based on atomic and ionic resonances, the acceptance 
bandwidth is lower than a few GHz [22], which definitely 
lies within the present assumption. 



3. Spatial dependance of the nonlinear susceptibility 

Although it is possible to calculate 70 m for an ar- 
bitrary R(z), this function is generally chosen so that 
one term 70^ is as high as possible and is the only op- 
timally phase-matched one. This is achieved by period- 
ically poling the nonlinear susceptibility with a period 
A. In this case, rh should be chosen such that L = mA, 
and = 2/tt, is the first term of the Fourier expan- 
sion of a square periodic function. In the following, we 
will consider the case of a collinear interaction in a pe- 
riodically poled crystal with an effective susceptibility 
Xeff = x(u; So ,a;£ )2/7r, corresponding to specific polar- 
izations of the pump, signal and idler beams. The case 
of the homogeneous crystal could also be described using 
an infinite period and R^i = 1. 



4- Quasi- degenerate down- conversion 



In the following, we will restrict the process to quasi- 
degenerate down-conversion, that is \Slj\ <C u Po with 
Suj = u Sn - uj in . 



The phase mismatch then reduces to 



AK(K s ,Ki) w Ak 



2k 



5uj \ n, 
uj P0 J n, 



p 



1+ A^npjK*!: (21) 



where Ako = k SQ + ki — k Po + ^ is the longitu- 
dinal phase mismatch. This together with the narrow 
bandwidth assumption allows for omitting the frequency 
dependence of the functions C^o, describing the target 
mode. It can therefore be replaced by the unique function 
0{) = W /2,o ~ O u ^_() ~ O u 



,0- 



D. Figures of merit 

In our general configuration of Fig. 1, only the spec- 
trally and spatially filtered pairs with a probability P2 
will give rise to measured coincidences between the two 
channels. It is also important to be aware of the to- 
tal brightness, accounting for all down- converted pairs in 
the spectral bandwidth under consideration Pq . The pair 
coupling efficiency is defined as the ratio T2 = P2/P0. 
This energetical figure of merit, however, cannot describe 
the loss of coherence due to spatial filtering, related to the 
ratio of coincidences to single counts. That is why we are 
interested in the probability r 2 | 1 to have the idler pho- 
ton coupled to the target spatial mode when the signal 
is coupled to that mode (or vice versa). It is sometimes 
called heralding ratio or conditional efficiency, measuring 
the ability of one photon to herald its twin for single pho- 
ton source applications. It is also very useful when the 
source is to be used as an entangled photon pair source 
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as a quality figure of merit. It will be evaluated using 
the single-photon coupling probability Pi of having at 
least one photon transmitted by the filters. In the fol- 
lowing subsections, a detailed calculation of each of these 
parameters is given. 



1. Source brightness 

The coincidental presence of both photons in the target 
spatial mode over a whole pulse duration is based on the 
two-photon state \^ 2 ) (Eq. (18)) and is given by 



(^2^2). 



(22) 



To make further calculations, we take the limit of an 
infinite quantization volume, in which sums over modes 
£ are replaced by integrals over the wavevectors k in the 
usual way [43, 44] with discrete variables like k^, uj£ and 
operators like o U£ changed into their continuous counter- 
parts k, cj(k) and o(o;(k)). Moreover, the narrow band- 
width assumption allows for calculating integrals over k 
independently of integrals over lj, since n mostly varies 
with the wavevector angle when uj is restricted to a small 
interval. Equation (22) gives, under the previous assump- 
tions 



p 2 = \cfn 2 K 2 , 



(23) 



where 



C = ie 



i(n' a u>e +n'iUi )zo/c ^pXeff^so^io ^\ 

y Seoc 3 n p n s ni 



is a constant with n' s (resp. n[) the refractive index in 
the medium where Go describes the target mode. 
The functions 



^2 ^ 

and 
K 2 



duji 
~2tT 



T p (uj s +cji -u Po ) 

2 

x F(uj s - (jj SQ )F((jji - uj io ) (25) 



[f d 2 n s 
J J (2^) 2 



d 2 t 



(27T) 



~2 Sp(K s -\- K-^, 0) 



x Oo(K a ,zo)Oo(Ki,zo)e ° { < kp ° < kp ° } 

x L sine 1 — 

V 2 



(26) 



describe respectively the spectral and spatial dependence 
of P2 . They must be both maximized in order to optimize 
the brightness. 

The function ^2 involves the coupling of the temporal 
dependence of the pump pulse with the spectral filters. It 
has the dimension of a frequency and can be interpreted 
as the effective source bandwidth. Using an adequate 



change of variables, it can be written as a convolution of 
the three spectrally-dependent functions and, as such, it 
is maximum when the filters are tuned so as to satisfy 
the energy conservation uj Sq + uji c 



U P0 



n 2 



|T P | 2 *|^| 2 *|^|V S 



u P0 ) (27) 



In the following, we will assume that this energy conser- 
vation condition is satisfied. 

The dimensionless function K 2 takes into account the 
spatial interferences caused both by phase matching and 
coupling to the target mode. Its maximization will re- 
quire a numerical optimization (see IIIC2). 

The factor \C\ 2 appears then as a spectral probability 
density. 



2. Pair coupling efficiency 

The pair coupling efficiency I^ is the ratio of the 
spectrally and spatially filtered pair probability (source 
brightness) to the spectrally filtered pair probability (to- 
tal brightness) 



P2 

Po' 



(28) 



The total brightness is calculated from the spectrally 
filtered, free space state |^o) defined by Eq. (16) 

Po = (^o|^o). (29) 

Using the same method as for P 2l Eq. (29) leads to 

P = \C\ 2 n 2 K , (30) 

where Kq is a dimensionless function describing the spa- 
tial dependence and defined by 



1 d 2 Ki 



Sp(ft s -\- 0) 



x Lsinc , 

V 2 

The pair coupling efficiency is then equal to 

r K2 



(31) 



(32) 



3. Heralding ratio and single-photon coupling 

Evaluating the heralding ratio requires calculating the 
single-photon coupling probability defined by 

Pi = (V>i|V>i> (33) 

for the state 

IV-i) = E ^s^i)l i T\Q4\Qol l a\ i \Q), (34) 
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where 



(35) 



corresponds to a spectral filtering of the signal and no 
filtering for the idler [45], and 



,ts Z0 



(36) 



describes the coupling of the signal photon only into the 
target single mode. 

In the same way as for P2 and Pq, 



Pi = ICffiitfi, 



where 
Oi = 



J 2ttJ 



duji 



d 2 Hi 



(27T) 5 



d 2 K, 

W) 



2 S p (n s + K,i,0)do(K s ,z ) 



_ i20 ^^ . AKL 
xe n s k Po Lsmc 



(37) 



(38) 



(39) 



As i^2, K\ requires numerical computation, whereas 
Vti reduces to 



/ 



dio 
2n 



\H« 



(40) 



which is proportional to the filter bandwidth itself. 

The heralding ratio i.e. the conditionnal probability 
to couple the second photon to the target mode when its 
twin has been coupled, is defined by 



r2|i — 



K2 



(41) 



It is useful to define also the single photon coupling 
efficiency 



r - Kl 



(42) 



These two coefficients are related to the pair coupling 
efficiency in the following way 



r 2|i Ti- 



(43) 



Note that r 2 7^ Tf, because of the spatial correlation 
between the two photons, in the same way as, because of 
the energy correlation of the two photons, the calculation 
of 0,2 involves a convolution rather than a product of the 
filtering functions of the signal and idler photons (see Eq. 
(27)). 

In order to get a physical insight into these results and 
show how this description can be used as a predictive tool 
for the design and optimization of entangled photon pair 
sources, the following section is devoted to a numerical 
calculation corresponding to the particular case of the 
experiment described in Section IV. 



III. NUMERICAL OPTIMIZATION OF A 
NARROW-BAND FIBERED SOURCE PUMPED 
BY A GAUSSIAN BEAM 



A. Gaussian pump beam, Gaussian target mode 



We assume the pump beam to be Gaussian with a waist 
radius wq and Rayleigh length zr 



S p (p,0) 



-\p?l< 



(44) 



(45) 



Using the fact that k Po /2 = zr/wq, the phase matching 
function becomes 



AK w Ah 







4z R 



2(1-^)^-2(1 + ^)^1 



Suj 



(46) 



The target mode acting as a spatial filter is defined by 
the profile of a Gaussian mode at its waist of radius ao 
located at z$. This can be the transverse mode of a fiber 
or its image by a lens collection system 



O (p,zo) = A / — o e -l^l 2 /«§ 



O (k,z ) = J2nal e"N a »/\ 



(47) 



(48) 



Nondimensionalization 



The coincidence probability P2 given by Eq. (23) can 
now be detailed. It is useful to separate the fixed param- 
eters from the configuration of the experiment that can 
be optimized 



gpXlfiLAuvu, 



8e c 4 n s rii 



O2 ,~ K 2 (t v 



k Po L 



(49) 

where Aujp is the filter bandwidth, and the terms 
Q2/AUJ-F and K2/(k Po L) depend only on the experimen- 
tal configuration described by the following dimensionless 
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parameters : 

4 In 2 

= — (relative pump bandwidth) 

L 



(normalized target mode waist size) 
(longitudinal collection offset) 










(1 r\ 


a = 






w 


c = 




L 


= 


Ak L 




w 


Vs = 


2 Ks 




w 


<Pi = 


2 Ki 



The parameter At p is the pump pulse duration. Eq. 
(49), that will be used to calculate the pair coupling 
efficiency and the heralding ratio, has a particular im- 
portance in itself, as it quantifies the absolute coupled 
brightness and its dependence on various experimental 
parameters. It will be discussed and compared to other 
reported work in Sec. V. The following section is dedi- 
cated to the theoretical determination of the experimen- 
tal configuration that maximizes this brightness. 




100 



FIG. 2. Effect of the relative pump linewidth 6 on the spectral 
transmission factor Q2 / Acjf 



for a photon and its twin to be transmitted is optimal. 
When the pump linewidth increases, some idler photons 
at uoi = u p — uj s and their corresponding signal photon at 
uj s are not frequency-symmetric with respect to the filter 
center frequency (lj p ^ u Po ), which is necessarily detri- 
mental. This effect becomes significant for Acu p > Acjf 
(5 > 1). Although reducing the pulse duration At p in 
order to increase the repetition rate does increase the ef- 
fective photon pair rate, it is desirable to compromise in 
order to keep ^2 close to its maximum value. 



C. Optimization of the source brightness 

1. Spectral optimization 

The optimisation of Eq. (25) must generally be made 
numerically, but the influence of the relative pump band- 
width S can be developped in a fully analytical way if the 
pump temporal envelope and the spectral filters are both 
Gaussian 



T P (t) 



41n2 



-21n2t 2 /Atp 



In 2 

-21n2w'' : /Aw 



e -^ 2 At=/(81n2) 



(51) 

(52) 
(53) 



where At p and Acjf are full widths at half maximum 
intensity. The coincidence probability then depends on S 
in the following way 



(S) 



/ 7r/(81n2) 
1 +T 



(54) 



This dependence of the spectral factor of the source 
brightness is plotted in Figure 2 as a function of the 
relative pump bandwidth. It shows that the maximum 
value is asymptotically reached when the pump beam 
is monochromatic. In this case, the joint probability 



2. Spatial optimization for degenerate down- conversion 

In the following, we restrict ourselves to the case of a 
frequency-degenerate down- conversion (u s = u^), which 
corresponds to the experimental setup described in sec- 
tion IV. The dimensionless term to be optimized has the 
following expression 



K 2 
k Po L 

with 



(f,a,C<Po) = -~£a 4 



d 2 fi Q2 



Q 2 = expj - \cp s + <^| 2 } 

xexp{-a 2 (|^| 2 + |^| 2 )} 

xex pi {-4,c(> s | 2 + |^| 2 )} 



(55) 



(56) 



x smc 



{f + i 



\<Ps + <Pi 



n s 



m J J 



Using polar coordinates, we can use the following map- 
ping : 



rOO pOO /»27T 

— > 2tt / p s dp s / pidpi / d(0 s - 6i) Q' 2 (57) 
Jo Jo Jo 
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FIG. 3. Optimization of the spatial filtering factor K2/ \k P0 L) 
with respect to the normalized target mode waist a and the 
longitudinal phase mismatch cpo for three values of the focus- 
ing parameter : £ = 0.1 (a) , £ = 1 (b) , £ = 10 (c) . 



FIG. 4. Maximal value of the spatial filtering factor 
K2/ (k PQ L) (a) and associated optimal values of the normal- 
ized target mode waist a (b) and the longitudinal phase mis- 
match (fo (c), for various values of the focusing parameter 



with 



h = ex P 

x 



x 



{-(l + a»)(p2 + p?)} 
exp I - 2p s pi cos(0 s - f ) j 

-"{- 4 «(2m + 3w)} 



x sine 



(1-2^ + (1-2^ 



2p s piCos(9 s -0i) j. 



(58) 



In this way, the quadruple integral turns into a triple 
integral which is numerically evaluated using an adapta- 
tive 3D quadrature algorithm [46]. 



Figure 3 represents, for three different focusing param- 
eters £, the value of the spatial filtering term K2/(k Po L) 
as a function of parameters a and ipo . It shows that there 
exists a unique couple (a opt , (/?o opt ) that allows reaching 
the maximum K% pt /(k Po L) for a given value of £. More- 
over, that maximum varies with £ and the accuracy of 
(a opt , (/?o opt ) is found to be more critical for low £. In- 
deed, when the focusing of the pump beam increases, 
phase matching can only be satisfied in an average way, 
because of the large range of emitted angles, and a opt 
results from a compromise between collecting weakly di- 
vergent photons with high efficiency and reducing this 
efficiency to collect more strongly divergent photons. Op- 
timizing the pump beam for photon pair collection indeed 
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consists in finding the focusing parameter for which these 
two compromises offer the best performance. 

Note that in Figure 3, the normalized longitudinal off- 
set of single-mode collection ( has been set to zero. We 
have checked that this choice gives optimal results. The 
optimality of ( = is due to the symmetry of the prob- 
lem in our particular choice of a Gaussian pump beam 
and Gaussian target spatial mode ; this might not be the 
case in other circumstances. 

In order to find the focusing parameter £ that gives 
rise to the highest collected brightness K% pt /(k Po L), we 
have performed optimizations similar to that of Fig. 3 
for values of £ ranging from 0.03 to 40, that is for waist 
radiuses from ~ 200 to ~ 5 /am for a red pump in common 
crystals, for instance. Knowing that £ = is optimal in 
the whole range, we have plotted on Fig. 4 the values 
K° pt /(k P0 L), a opt and ^ O opt as a function of f . 

According to Fig. 4b, the optimum normalized target 
mode waist a opt does not vary very much over this large 
range of focusing parameters. Starting with an approx- 
imate matching of the target mode to the pump waists 
a o P t ~ i a t \ ow focusing, our calculations exhibit a slow 
increase, showing that when focusing gets stronger, it is 
preferable to collect the weakly divergent photon pairs, 
since modes of larger waist sizes have smaller numerical 
apertures. On Fig. 4c, the optimum longitudinal phase 
mismatch (po° pt increases with the focusing parameter, 
as to compensate for the transverse mismatch caused by 
the strong focusing. 

Fig. 4a shows that the value of the focusing parame- 
ter giving the highest value of K^ pt /(k Po L) is £ = 2.84, 
for which the longituninal phase mismatch is ipo = 3.2. 
These values correspond to the Boyd and Kleinman con- 
ditions [47] for a highest second harmonic generation ef- 
ficiency with Gaussian beams. Indeed, when evaluating 
the brightness produced in a single Gaussian mode, the 
collinear degenerate down-conversion is a process sym- 
metric to second harmonic generation. 

However, down- converted photons are not entirely pro- 
duced in a single mode. That is why it is important to in- 
vestigate the optimization of the pair coupling efficiency 
r 2 defined in Eq. (28). 



D. Optimization of the pair coupling efficiency 

The total source brightness must also be calculated, 
using the dimensionless variables defined previously : 

Po = ^^ .^(i).^^), (59) 



Se c 4 n s ni 



k Po L 



where 



K = 2_ 

k po L 7T 3 ' 



(60) 
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FIG. 5. Pair coupling efficiency T 2 (a) compared to the spatial 
factor Ko/(k PO L) of the pair production probability (b), for 
various values of the focusing parameter £. Ko/(k PO L) is eval- 
uated for the longitunal phase mismatch ipo that maximizes 

r 2 . 



with 



Qo = expj - \<p s + <^| 2 } 



x smc 



\<Ps + <Pi 



^\<Ps\ 2 



1> 



(61) 

I 2 ]} 



which can also be reduced to a 3D-integral. 

The pair coupling efficiency T2 (Eq. (32)) is calculated 
using Eqs. (55) and (60). It is plotted in Figure 5a, 
along with the dimensionless Ko/(k Po L) proportional to 
the total brightness Po (Fig. 5b), as a function of the 
focusing parameter £. The total brightness does not vary 
much for £ > 1, therefore the coupling efficiency follows 
the same tendency as the single-mode brightness (Fig. 
4a). However, the precise optimal value of £ is no longer 
£ = 2.84 but is closer to 2. For a red pump in common 
crystals, this means a waist radius difference of 20 %. 
Under different assumptions from the ones considered in 
this section, however, the difference might be greater. 
When power efficiency is critical, as for sources intended 
to be compact, it can be useful to have in mind that 
to collect most of the generated photons, optimizing on 
the mere source brightness is not optimal with respect to 
pump power consumption. 
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E. Optimization of the heralding ratio 

The heralding ratio r 2 |i (Eq. (41)) requires the com- 
putation of the single-photon coupling probability Pi , (cf. 
Eq. (37)) : 



8eoc 4 n s iii 



Ki 

Auj-p K 7 k Po L v 



(62) 



where only depends on the shape of the filter itself 
and where 



with 

Qi =exp{ - IVs + V'il 2 } 
xexp{ -a 2 |<y? s | 2 } 



(63) 



(64) 



X 



exp i{-4eC^|^| 2 } 
mc < - 



x smc 



2^|^| 2 -2 



/7 ; J J 



which can be reduced to a 3D-integral as for Q2 and Qq. 

The heralding ratio r 2 |i = K 2 /Ki is plotted in Fig. 
6 as a function of parameters a and ip$ for three values 
of the focusing parameter £ = 0.1, 1, 10. Contrary to the 
brightness depicted in Fig 3, the heralding ratio can reach 
a value close to 1 on the whole range of £. A large range 
of (fo is compatible with this maximum, but the overlap 
with the range leading to a high brightness is small. On 
the contrary, the tolerance on a is relatively low. When 
the results of P 2 and r 2 | 1 are both taken into account, the 
theory gives indeed useful information about the target 
mode waist for which the collection should be optimized 
with respect to ip$ (through the crystal temperature), to 
find the best compromise between the brightness and the 
heralding ratio. 

The details of an experiment that enabled a validation 
of our model is given in the following section. 



IV. EXPERIMENTAL VALIDATION 

A. Experimental setup 

The experimental setup used to validate the theory 
is depicted in Fig. 7. SPDC is generated by focusing 
a pulsed pump beam at 782 nm in a periodically poled 
lithium niobate (PPLN) crystal of length L = 2 cm with 
a poling period A = 19.34/im. The mean pump power 
is 5 mW and the 25 ns Gaussian pulses (FWHM) are 
Fourier transform limited with a 2 MHz repetition rate. 
The spatial profile of the pump beam is also Gaussian. 
On- axis fluorescence around 1564 nm is collected into a 
telecom optical fiber through the lenses L c (achromatic 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 

a 

FIG. 6. Optimization of the heralding ratio r 2 |i with respect 
to the normalized target mode waist a and the longitudinal 
phase mismatch (po for three values of the focusing parameter : 
£ = 0.1 (a) ,£ = 1 (b),£ = 10 (c). 



doublet) and Li (asphere). The same low bandwidth fil- 
ter ( Acjp = 27r x 75 GHz) is used for both signal and idler 
photons so that the source is operated at the degeneracy 
frequency (uj Sq = c^ = uj Po /2). A balanced coupler is 
used to split photon pairs with 50 % efficiency and the 
photons are detected on paths A and B. 



B. Experimental method 

In order to explain how the heralding ratio r 2 |i can 
be determined experimentally, let us show its relation to 
measured parameters : 



r 2|l 



K2 



(65) 
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Pump 




FIG. 7. Experimental setup : A pulsed pump laser at wave- 
length 782 nm is focused in a periodically-poled lithium nio- 
bate crystal (PPLN) with a lens L p of focal length f p . Down- 
converted photons at 1564 nm are coupled into an single mode 
fiber through an optical system composed of lenses L c and Li 
of respective focal lengths f c and fi. Spectral filtering is per- 
formed via a DWDM add-drop filter of bandwidth 75 GHz, 
and photons are split with 50 % efficiency towards detectors 
A and B using a balanced fibered coupler. 



The experimental parameter Pab corresponding to the 
calculated P 2 is the measured coincidence probability per 
pulse from which accidental and noise coincidences are 
substracted, and Pi is related to the measured counts on 
channel I from which dark counts have been substracted : 

Pab = T A T B P2 Pi = 7}Pi (66) 

where 7} is the overall tranmission of channel I. We 
obtain : 



F2|l = 



Hi Pab Ti 
n 2 TaTb Pi 



(67) 



The heralding ratio r 2 | 1 can hence be determined from 
the measurements of counts and coincidences, provided 
the insertion losses have been previously determined [48] . 
It is then possible to validate its dependence with respect 
to the pump focusing parameter £ and the normalized 
target mode waist a. 

The variation of £ was obtained by changing the lens 
L p focusing the pump beam into the PPLN crystal. For 
each value of the focal length / p , the waist of the pump 
beam was measured, allowing the determination of zr 
andf = L/(2*fl). 

Counts and coincidences were then measured using 
various focal lengths fi of the lens focusing the SPDC 
beam into the fiber. The value ao of the image of the 
fiber waist in the crystal was determined using the mag- 
nification factor fd fi of the collection system (composed 
of the lenses L c and Li) in order to determine a = ao/w^. 
Let us note that each data point requires changing the 
focusing lens, realigning the setup, and successively opti- 
mizing the collection with at least five different injection 
lenses, with the phase mismatch (crystal temperature) as 
an additionnal degree of freedom. 



C. Experimental results 

The experimentally determined r 2 | 1 can then be com- 
pared to its theoretical value when the normalized target 
mode waist a is varied. Figure 8 shows two examples : 
the pump focusing is kept constant with f p = 100 mm 
(f = 0.76, Fig. 8a) or f p = 50 mm (f = 2.7, Fig. 8b). 
r 2 |i is normalized and plotted as a function of a for ( = 
and the adequate value of the remaining unknown dimen- 
sionless parameter cpo is found by horizontally fitting the 
theoretical curve to the experimental data. For £ = 0.76, 
(fo is found to be around 2.0 while it is around 3.2 for 
£ = 2.7. These values of cpo are also in agreement with 
the theoretical predictions corresponding to the bright- 
ness optimization. Indeed, although each plotted val- 
ues of r 2 |i was obtained after optimizing this figure of 
merit itself, the experimental starting point was a pre- 
liminary optimization around a maximum source bright- 
ness. When r 2 |i and P 2 have a common optimum in 
the (a, (fo) space, it is not suprising to converge close to 
(a opt , (po opt ) (optimum of P 2 as shown in Fig. 4 as a 
function of £) when optimizing with respect to r 2 |i- 

The very good agreement of our theory with the ex- 
perimental results confirms its validity. For practical ap- 
plications it is remarkable that only the measurement of 
the pump waist is required for the optimization of the 
SPDC source, the choice of the collection magnification 
being directly derived from the calculations presented in 
our analysis, which leaves the remaining parameter cpo to 
a heuristic, using the mere crystal temperature. 

The absolute value of the heralding ratio was also de- 
rived from the preliminary measurements of Ua = 0.026 
and Ub — 0.024. In Fig. 9a, the experimental absolute 
value of the heralding ratio r 2 M is plotted as a function of 
the focusing parameter £. The result is almost constant, 
as predicted by the theory, but around 30 % below the 
expected optimal value that is close to 100 %. This dif- 
ference is probably due to imperfections in the Gaussian 
pump beam and aberration in the optical system used to 
eliminate pump photons and to collect down- converted 
photon pairs into the fiber. 

Let us also remark that the maximum measured value 
of r 2 |! =0.7, although not optimal, is to our knowledge 
impossible to reach by coupling PPLN waveguides into 
single- mode fibers, due to the mode mismatch between 
the rectangular profile of the guided mode and that - 
circular- of a fiber. This makes focused SPDC in bulk 
crystals probably more suitable than waveguide crystals 
when a high state fidelity is critical and limited by other 
constraints. For instance, the complex spectral filtering 
required for quantum memories associated with its lim- 
ited storage-retrieval efficiency makes the coupling effi- 
ciency a parameter that is more than ever desirable to 
maximize. 

In summary, we have confirmed 1) the dependence of 
the heralding ratio on the target mode size, 2) its quasi- 
independence on the focusing parameter, and 3) that our 
theory predicts with a very good precision the optimal 
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FIG. 8. Comparison of theoretically (solid lines) vs. experi- 
mentally (dots) determined heralding ratio r 2 |i, as a function 
of the normalized target mode waist a for two different focus- 
ing parameters : £ = 0.76 (a) and £ = 2.7 (b). 



target mode size for a particular focusing parameter. 



V. DISCUSSION 
A. Comparison with other works 

As mentioned in the introduction, our work covers 
issues that were addressed by other authors with vari- 
ous assumptions or methods. The two figures of merit 
that can provide a relevant comparison are the absolute 
brightness and the coupling efficiency. We will also dis- 
cuss a few reported experimental validations. 

About the absolute brightness, Ling et al. [49] were 
among the first to propose an expression of the absolute 
photon pair rate collected into a single Gaussian spa- 
tial mode. Apart from the fact that they consider the 
case of a continuous pump ((5 = 0), our calculations are 
consistent dimension-wise. Moreover, everything else be- 
ing equal, a spatial filtering factor can be identified in 
their analysis: its dependence with the focusing param- 
eter and the longitudinal phase mismatch is of the form 
£sinc 2 ((/?o/2), which is indeed the result we find for the 
specific case of low focusing < 1) and a = V% cor- 
responding to Ling's assumptions : thin crystal, negli- 
gible pump diffraction, equal Rayleigh lengths for pump 
and down-converted photons. However, under our less 



FIG. 9. (a) Measured heralding ratio F 2 \i plotted for 4 focal 
lengths f p = 150, 100, 75, 50 mm, as a function of the focus- 
ing parameter £ = 0.44,0.76,1.05,2.70 . (b) Corresponding 
experimental (circles) and calculated (solid line) value of the 
normalized target mode waist a for each value of r 2 |i 



restrictive assumptions, we have shown that a = \[2 is 
optimal for £ = 2.84, whereas for low focusing, a < 1. 

More recently, the absolute value of the source bright- 
ness was calculated by Mitchell [37] in the Heisenberg 
picture, under assumptions similar to ours. This inves- 
tigation has been restricted to a monochromatic, contin- 
uous pump field, and a direct expansion of the Hamil- 
tonian onto the target Gaussian modes was performed, 
as opposed to our choice of investigating the single-mode 
coupling of free-space expanded fluorescence in the in- 
teraction picture. Our theory, which is developed in the 
Schrodinger picture, is nevertheless fully consistent with 
the results of Mitchell and we suggest in the next subsec- 
tion how to generalize it using the same framework. 

As for the coupling efficiency issue, let us recall that 
we investigated two parameters : the pair collection 
efficiency and the heralding ratio. Our results con- 
cerning the pair collection efficiency are consistent with 
that of Bennink [38], who found a linear dependence 
on the crystal length when the source bandwidth is 
much smaller than the phase-matching bandwidth, as 
opposed to Ljunggren et al. [35] who conclude on a 
V^-dependence. Let us point out that, as Bennink [38], 
we have taken into account the diffracting nature of the 
pump beam, and considered the longitudinal phase mis- 
match as a degree of freedom for optimization. As a 
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result, our optimization of T2 for Gaussian beams at 
the degenerate frequency is optimal close to the well- 
known Boyd & Kleinman conditions [47] : £ = 2.84, 
cpo = 3.2 and equal Rayleigh length for the pump and 
down- converted photons (a = y/2). An advantage of our 
framework is that the particular pump and target modes, 
filter shape, pump linewidth, etc. are only used in the 
very last step of the optimization, i. e. the computation of 
a mutliple integral. Up to then, our framework remains 
very general. 

As far as the heralding ratio is concerned, as opposed 
to the coupled brightness, it is found to be close to 1 
whatever the focusing strength, provided the target mode 
waist a and phase mismatch ipo are adjusted according 
to our theory. Interestingly, a subset of (a,(po) which 
maximizes the heralding ratio is generally close to the op- 
timal brightness, giving the configuration for an optimal 
general source performance. This conclusion is different 
from that of Benninck [38] who finds that a strong reduc- 
tion of brightness is necessary to achieve a high heralding 
ratio, highlighting what could be a fundamental trade- 
off for single- mode applications. This apparent disagree- 
ment could be due to our specific narrow-band assump- 
tion. Our experiments confirm that the narrow-band be- 
haviour follows our predictions, which is a positive result 
as far as the perspectives of SPDC for narrow-band ap- 
plications are concerned. 

As remarked by Bennink, very few among the reported 
theoretical works have been experimentally validated, 
probably because of the complexity of such experiments 
that involve multiple parameters to optimize simultane- 
ously and to measure precisely. Reported experiments 
were realized for a single set of non-optimal parameters 
[35] or in a configuration which is not in the scope of 
our theory [34]. Thanks to a careful reduction of sys- 
tematic errors, our experimental results show that our 
model could be used in order to optimize the experimen- 
tal source configuration. 



B. Extension to other source designs 

Although the calculation done in Sec. Ill uses a bulk 
periodically-poled crystal pumped with a single-mode 
diffractive gaussian pump beam, as required for success- 
ful comparison with our experiment, the general theoret- 
ical framework of Sec. II can be applied to other source 
designs. A few examples are given her after. 

First, it is straightforward to correct the results for 
non-degenerate down conversion, provided the frequency 
difference is small compared to the pump frequency. 

It is also possible to use a non-gaussian pump or target 
mode profile, by changing the functions S p and Oq re- 
spectively. If required, a second dimensionless parameter 
can be introduced apart from a in order to characterize 
the mode ellipticity. In this case, the 4D- to 3D-integral 
reduction used in Eq. (57) is impossible due to the ab- 
sence of cylindrical symmetry, making the computation 



longer. Let us note that for such high-dimensions inte- 
grals, a Monte-Carlo integration method could give faster 
results. 

As far as the crystal type is concerned, as mentionned 
at the end of Sec. II A, a non-periodically or multi- 
periodically poled crystal can be used by coherently sum- 
ming the first N components of the Fourier series ex- 

(2) 

pansion of the non-linear susceptibility % (z). When 

—(2) 

X (z) is not of high complexity, the expansion can be 
truncated to a low N, making the calculation not much 
more time-consuming. 

SPDC sources using a periodically-poled crystal with 
integrated waveguide are also commonly used. The cou- 
pling of photons down- converted in such devices could be 
modelled using our framework, by slicing the crystal into 
N sub-crystals, for each of which the state coupled to the 
waveguide mode is calculated. For a slice j = 1 • • • N cen- 
tered on Zj , the target mode is assumed to have its waist 
located at Zj . If diffraction is taken into account as in the 
case of bulk crystals, the pump profile S p (p, Zj) differs in 
each slice according to Eq. (10). However, the pump 
field is generally also guided in the integrated waveguide 
and the pump profile should be chosen constant, which 
is done by removing the diffraction term in Eq. (10). 
The final state for which figures of merit like brightness 
or heralding ratio are evaluated results from coherently 
summing the states originating from each slice, taking 
into account their respective phases. 

Let us remark that, in integrated waveguides designed 
to be single- mode for down- converted photons, the pump 
field is generally spatially multimode. This can also be 
taken into account by replacing the single-mode pump 
function S p by a mode expansion X^jLi^Pj- As f° r a 
general nonlinear susceptibility distribution, this only in- 
creases the computation time by a factor TV, which can 
be chosen according to the desired precision. 



VI. CONCLUSION 

In this paper, we have calculated the state of a photon 
pair produced by a narrow-band spontaneous paramet- 
ric down- conversion (SPDC) source with arbitrary pump 
spatial and temporal profile, and arbitrary filtering con- 
figuration. 

When applied to Gaussian modes, our theory is consis- 
tent with the most recently reported work where realistic 
assumptions have been made. Within our assumptions, 
no incompatibility is observed between a high brightness 
and a high heralding ratio. We believe this result could 
increase the interest of SPDC sources for narrow-band 
quantum information applications. 

Once validated under the Gaussian assumptions, our 
theoretical framework, which allows an extensive study of 
the source through many degrees of freedom, should allow 
the prediction of the best performance for more original 
source designs, including non-Gaussian pump and target 
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modes, guided or diffractive, single or multi-mode, with 
arbitrary nonlinear susceptibility longitudinal distribu- 
tion. 
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Appendix A: Conventions for Fourier transforms 
and convolutions 

All along this paper, we use the following Fourier trans- 
form conventions : 

/M = jdt e^f(t) <-> f(t) =j^- e-^/M 
/>) = |d 3 re- ik - r /(r)^/(r) = J-^ e ik r />) 
/ is normalized if : 

jdt\f{t)\ 2 = /^|/H| 2 = 1 

/ d3r|/(r)|2 = /(0 l/(k)|2 = 1 
Convolutions are defined as follows : 

f*g(t) = Jdt'f(t')g(t-t') 

f*g(r) = Jd 3 r'f(r')g(r-r') 

f*g(k) = J j^f(k')g(k-k>) 

Appendix B: Spatial and spectral filtering 

In this appendix we describe the effect of filters on a 
single photon state. Generalization to the case of photon 
pairs is made in Section II B. 

Although the action of filters will be described in the 
Fourier space, they are located at a particular position 
of the setup and a specific time lag [to j ^1] can be defined 
such that when t < to? the photon state is completely 
unaffected by the filter yet and when t > t\ , the effect of 
the filter on the photon state has been completed. After 
having introduced some tools in the first two sections, we 
will successively examine how spectral and spatial filters 
change a state |^(to)) m ^o a state \^{t\)). 



1. "Localized" photonic states and 
pseudo-wavefunction 

In this paper, we sometimes use localized photonic 
states introduced by L. Mandel [50]. Based on plane- 
wave states | lk) in a quantization volume V, a state de- 
scribing a photon localized around r is defined by : 

l 1 r) = E £ 7^l 1 k) ( B1 ) 
k * 

Then, a general photon state = J^k^kllk) can be 
written as 

M= j d 3 r^(r)|l r > (B2) 
Jv 

where -0(r) is a spatial pseudo-wavefunction for which if;^ 
are the coefficients of its Fourier series expansion in the 
quantization volume V. Such a wavefunction is only valid 
provided the volume V in which the localization probabil- 
ity -P(V), defined as follows, is evaluated is large enough 
(each dimension much larger that the wavelength) : 

P(V)= [ d 3 r\(l r M 2 =[ d 3 r\^(r)\ 2 (B3) 
Jv Jv 

This probability is naturally equal to unity in the quan- 
tization volume : 

P(V) = / d 3 r |V(r)| 2 = <V#> = £ |^k| 2 = 1. (B4) 



2. Spatial filtering 

a. Principle 

A spatial filter located on the propagation axis at 
z = zo is modelled as the coupling in the plane z = zq of 
the down-converted field into a single spatial mode de- 
fined by a function Oo;,o( r )- The index uo indicates that 
this spatial mode can be frequency-dependent for a given 
filter. 

Note that zo may be in [— L/2,L/2]. As an example, 
if the spatial filtering is done via an optical fiber, the 
effective location of the filter is where the lens collection 
system images the entrance of the fiber. 

To evaluate its transmitted component, the down- 
converted field initially described as a superposition of 
plane waves is better expanded on a particular set of or- 
thogonal modes {0 Wj j(r)}, one of which (Oo;,o( r )) being 
the spatial mode selected by the considered filter. 

If coupled to other modes, photons are supposed to be 
lost. As for spectral filtering, this leads to a mixed state. 

For instance, the fundamental mode can be the Gaus- 
sian mode of the Laguerre-Gauss basis, suited to single- 
mode fibers. Then j=o(r) describes a Gaussian beam 
of waist size equal to the field radius of a fiber mode. 
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If the image of the fiber entrance is located at z = zo, 
Oujj = o(p, zo) describes the Gaussian beam waist profile 
corresponding to the fiber transverse mode. 

A projection at z = zo into that single mode selects 
the photon state component which is transmitted by the 
filter. Once transmitted, the state is said to be in the 
spatial mode characterized by the annihilation operator 
o u . Let us note that whether the spatial mode for z > zo 
continues to be that of a Gaussian beam (as for a cavity) 
or is actually described by a guided propagation with a 
constant transverse profile Ocu,o(p, ^o) and propagation 
constant j3 u = nf(uo)uo/c like in a single mode fiber of 
effective refractive index n' [51] will not change the de- 
scription: 



O UJ j=o(p : z > z ) = w> o(p,2o)e : 



iPw(z—z ) 



(B5) 



In a quantization volume V = S x £, a field in the 
mode O u j(r) will be said to be in the following quantum 
state : 



\O uJ ) = -j=J<PrO uJ (r)\l I 



(B6) 



so that the states are normalized : (O ^ j\0 ^ jf) = 
Sjjf. Functions Wj j(r) are only transversally nor- 
malized : J d 2 pO^ j,(p, z)O UJ j(p, z) = 5jj> and 
J v d 3 rO: tJ ,(r)O 0jJ (r) = CS jr . 

b. Calculation of the transmitted field component 

In a quantization volume containing a crystal before 
z = zq and a spatial filter at zq, a one-photon field in 
mode I is in the state : 

|l,) = |d 3 r/(r)|l r ) (B7) 
+ e(z-z )Y^9 j (t)^f=ft)\U) (B8) 



where Q(z) is Heaviside's function and gj(f) are deter- 
mined by the boundary conditions at z = Zq : 



(B9) 



Remark that the wavevector takes into account the 
change of medium according to Snell-Descartes' law of 
refraction : 

k' _ ^in^ = i K _ Ik'lcos^ ^n> 
k |k|sin# k z |k|cos6> n 

where n' is the refractive index of the medium in which 
O 0Jij j(r) describes the spatial mode and n the refractive 
index in the crystal. 

Using the orthonormality of functions {O w j(r)}, one 
gets : 



9j {i)= i d 2 po: e A P ,z ) 



Ant- p f>ik' z £ -z 



Vv 



(Bll) 



where U£ j (k,, zq) is the Fourier transform of O^^j (p, zq) 
and * designates a complex conjugate. 

Starting from an initial state |\l/(to)) = lo- 
calized at z < zo, the state becomes for z > zq : 

l*(*l)> = J2J2^ eik ' zXZ0 ^ d ueA"t> Z °)\°^j) 
j t V6 

(B12) 

Ignoring components of the state which are not within the 
mode | O u£j o) transmitted by the filter, the transmitted 
one-photon state becomes : 

i*(*i)> = Ew^ /z °^;,o(«^)i^,o) 

+ ignored zero-photon states (B13) 

Note that if one of the functions {(D^-fr)} describes 
the spatial eigenmode of any other filter (e.g. a rectan- 
gular waveguide), the result is still valid. 

As for spectral filtering, extrapolation to two-photon 
states is straightforward. 
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